Does Pulsar B1757-24 Have a Fallback Disk? 
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ABSTRACT 

Radio pulsars are thought to spin-down primarily due to torque from 
magnetic dipole radiation (MDR) emitted by the time-varying stellar magnetic 
field as the star rotates. This assumption yields a 'characteristic age' for a 
pulsar which has generally been assumed to be comparable to the actual age. 
Recent observational limits on the proper motion of pulsar B1757~24, however, 
revealed that the actual age (> 39 kyr) of this pulsar is much greater than 
its MDR characteristic age (16 kyr) - calling into question the assumption of 
pure MDR spin-down for this and other pulsars. To explore the possible cause 
of this discrepancy, we consider a scenario in which the pulsar acquired an 
accretion disk from supernova ejecta, and the subsequent spin-down occurred 
under the combined action of MDR and accretion torques. A simphfied model 
of the accretion torque involving a constant mass inflow rate at the pulsar 
magnetosphere can explain the age and period derivative of the pulsar for 
reasonable values of the pulsar magnetic field and infiow rate. We discuss 
testable predictions of this model. 

Subject headings: Stars: neutron — pulsars: individual (PSR B1757-24) — ISM: 
supernova remnants 
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1. Introduction 



Isolated pulsars are spinning neutron stars whose observed spin rates gradually 
decrease with time. The age r of a pulsar is usually assumed to be equal to the timing (or 
characteristic) age tmdr derived by assuming pure magnetic dipole spin-down in vacuo, 
and the age is then given by tmdr = ~Q/{2Q) (e.g. [Manchester fc Taylor 19771) , where 
Q = 27c/P and Q = —2ttP/P'^ are the angular spin frequency and angular frequency 
derivative for a pulsar spin period P and period derivative P. Under this same assumption 
of MDR spin-down, the magnetic field strength of the pulsar is given by the formula 
( [Manchester fc Taylor 19771 ) 

5 = 3.2 X 10^^(PP)^^^G, (1) 



which is often assumed to be equal to the true field strength of the isolated pulsar (e.g. 
I'aylor, Manchester, fc Lyne 199^ ). 



PSR Bl 757-24 is a 0.125 s radio pulsar which appears to be associated with the 
supernova remnant G5.4-1.2 ( jUaswell et al. 1987[ ). The pulsar is surrounded by a 
compact radio nebula having a cometary morphology with a tail extending back into the 
supernova remnant ( [I'rail, Kassim, fc Weiler 1994[ ) - strongly suggesting that the pulsar 
was formed in the supernova which produced G5.4-1.2 ([Manchester et al. 1991 ). Given the 



temporal parameters of PSR B1757-24 (P = 0.125 s and P = 1.28 x 10-^^ s s^^ [laylor 



[Manchester, fc Lyne 19"93[) , tmd_r, = 16 kyr and P = 4xlO^^G for the pulsar. Assuming 



that Tmdr is similar to the pulsar's its true age, the transverse velocity implied by the 
pulsar's displacement from the apparent center of the G5.4-1.2 is greater than 1500 km 
s^i ( [Gaensler fc Frail 200q ; [Frail fc Kulkarni 1991[) . Observations of PSR Bl 757-24 taken 



six years apart, however, failed to detect the expected proper motion from the pulsar, 
yielding a distance-independent lower limit on the age of PSR B1757-24/G5.4-1.2 of 39 kyr 
( paensler fc Frail 2000|) . This is more than a factor of two greater than the pulsar's MDR 
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characteristic age. 

The discrepancy between the proper motion age and the MDR characteristic age of 
PSR Bl 757-24 suggests that the spin-down of the pulsar is not due purely to MDR but 
also has significant contributions from other sources of torque. Istomin (1994) considered 
a model for PSR B1757-24 in which the pulsar was interacting with dense plasma in the 
shell of G5.4-1.2, causing an increase in the torque at the light cylinder. Here we consider 
another possible source of extra torque on the neutron star - from a disk of material 
accreted from ejecta produced in the supernova explosion. These fallback disks may be 
roughly divided into two categories: "prompt" and "delayed" . Prompt disks may be formed 
from ~ 0.001 — O.IM© ( [Michel 1988 ; Lin, Woosley, fc Bodenheimer 19OT ) of ejecta material 



soon after the initial core collapse in a type II supernova explosion ([Woosley fc Weaver 



|1995|) . Formation of such prompt disks is probably limited to < 7 days after the core 
collapse because of heating of the ejecta by ^^Ni decays ( [Chevalier 1989| ). Delayed disks may 
form years after the explosion from ejecta decelerated by radiative cooling ([Fryer, Colgate, 



Pinto 1999| ) or by a strong reverse shock ( [Truelove fc McKee 1999[) caused by the primary 



supernova blast wave impinging on dense circumstellar material from the pre-supernova 
stellar wind ( [Gaensler 1999[ ). Whether or not a neutron star accretion disk can form shortly 
after a supernova explosion depends on the opposing forces of the pulsar MDR wind and 
the pressure of the hot and turbulent environment shortly after the explosion. Since the 
latter is highly uncertain ( [Woosley fc Weaver 1995[ ), for the purposes of this Letter we 



assume that a disk can form around a neutron star under these conditions and explore the 
implications for PSR Bl 757-24. 



2. Spin-down From Accretion Torques 
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An accretion disk around a magnetized neutron star can exert a spin-down torque on 
the star if the mass inflow rate is low and the magnetic field is strong. Quantitatively, 
this condition is met when the Keplerian co-rotation radius Rc = 1.7 x lO^P'^^^ cm is less 
than the magnetospheric radius Rm = 4.6 x 10^Bi2^mi^^'^ cm, where m = lO^^rhis g s~^ 
is the mass inflow rate onto the magnetosphere and B = 10^^i?i2 G is the strength of the 
neutron star dipole magnetic fleld (see e.g jl'rank. King, fc Kaine 1992[ ). Here and elsewhere 



we assume a neutron star mass and radius of M* = IAMq and -R* = 10 km, respectively. 
When Rm > Rc, infalling material is stopped at the magnetosphere by a centrifugal barrier 
which prevents accretion onto the neutron star surface. Instead, the infalling material 
may be accelerated away in a wind which carries away angular momentum from the 
magnetosphere and hence the neutron star itself. This "propeller effect" ( piarionov fc 



Sunyaev 1975| ) spin-down mechanism has been invoked to explain the behavior of some 



Galactic accretion-powered x-ray pulsars ( |Cui 199^ ; [Zhang fc Yu 1998|) , the spin-evolution 



of anomalous x-ray pulsars (AXPs: [van Paradiis, Taam, fc van den Heuvel 1995 ; Chatter] ee 



Hernquist, fc INarayan 2(]UD| ) and soft gamma-ray repeaters (Alpar 2(J(JC; Marsden et al. 



|2000| ), and the optical and infrared spectra of some radio pulsars ( [Perna, Hernquist, fc 
Narayan 2000| ). 



The pulsar B1757-24 is in the propeller regime (-Rm > Rc) for mass inflow rates of 
m < 5 X 10^'' g s~^ and a canonical ( [Manchester fc 'lay lor 19T7| ; [laylor, Manchester, fc 



Lyne 19931) neutron star magnetic fleld of ~ 10^^ G. Since the characteristic age for an old 



[> 10 kyr) pulsar depends only weakly on the initial spin period ( [Manchester fc 'laylor 



1977|) , we neglect the different formation times of the delayed and prompt fallback disks. 



and assume that the spin-down evolution of the pulsar at all times is determined by the 
combined torque from both the MDR and the accretion disk. The spin-down rate due to 
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radiation from a rotating magnetic dipole is ( |Manciiester fc 'lay lor 19T?|) 

= :pr^ — , (2) 



where = ^M^,Rl is the neutron star moment of inertia. The spin-down rate due to the 
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propeller torque vanishes gradually (IChatterjee, Hernquist, fc Narayan 20001 ; |Alpar 2000| ) 



as the star approaches spin- equilibrium (Rm = Rc) at a spin period P,, = A.7B%' 
The spin-down rate due to the propeller torque is simply ([Menou et al. 1999| ) 

si. = (i - 1-) , (3) 



where A; is a positive constant of order unity ( |Wang & Robertson 1985[ ) and il^q = 2ti / Pf.q. 



Assuming a constant mass inflow rate and dipole magnetic field, the timing age for a final 
spin period P = 27r/fi is given by 

'^comb / rS rS ' ^ ' 

where VLq is the initial angular frequency. A more realistic expression for the propeller 
torque would incorporate a time dependent mass inflow rate ( |Cannizzo, Lee, fc Goodman 



1990 ) in Equation (^, as m should decrease in time as the disk dissipates. For the model to 
be correct, however, a disk must still be present around the pulsar, because otherwise the 
MDR timing age would be greater than the true age. Therefore the m used here may be 
thought of as a time-averaged value of the mass inflow rate. In addition, the effect of the 
propeller flow on the MDR torque ( [Roberts fc Sturrock 19731) is not taken into account. We 



plan on incorporating both of these effects in future work. 

A contour plot of the pulsar B1757-24 timing age Tcomb for various values of the 
magnetic field strength B and mass inflow rate m is show in Figure 1. The characteristic 
ages were calculated using Equations and assuming P = 0.125 s, Pq = 10 ms, and 

k = 1. In this simple model, the allowed values of B and m for pulsar B1757-24 lie on the 
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heavy solid line corresponding to the observed P and P. In addition, the shaded regions 
of parameter space in Figure 1 are excluded by the lower limit on the age (right shaded 
region: paensler fc Frail 20001) , and the necessary condition Rm > Rc (shaded region in 
upper left corner). From Figure 1, we find that values of 2 x 10^^ < B < 1.4 x 10^^ G, 
7 X 10^^ < m < 7 X 10^'' g s^^, and 39 < Tcomb < 60 kyr are consistent with the lower limit 
on the true age of 39 kyr ( [Gaensler fc 1^'rail 2000| ) and the present-day spin-down rate of the 
pulsar. 



3. Discussion 



In the context of this model, the required mass inflow rate for PSR Bl 757-24 overlaps 
the range of m inferred from accretion-powered neutron star systems ([Bildsten et al. 1997| ). 



Radio pulsations are not seen from accretion-powered x-ray pulsars in binary systems 
( Pender et al. 1996|) , which implies that the emission mechanism responsible for radio 
pulsations may be quenched by matter near the polar caps of the pulsar. This is not 
a problem for the PSR B1757-24 model, however, because in propeller sources most of 
the matter is ejected before it has a chance to reach the polar cap and quench the radio 
emission. Radio emission may be suppressed for propeller sources closer to spin equilibrium 
than PSR Bl 757-24 (e.g. the AXPs: phatterjee, Hernquist, fc Narayan 2000|) , however, 
because as equilibrium is approached more matter will be allowed to fall onto the neutron 
star surface. The total mass required to fuel the propeller spin-down over the lifetime (so 
far) of pulsar B1757-24 would be mr ~ 10^® — 1O~^M0 - a tiny fraction of the total amount 
of ejecta in a typical Type II supernova explosion ([Woosley fc Weaver 1995| ). 

This hypothesis can be tested by multiwavelength observations. In this model, 
the total emission from the pulsar would be due to the propeller wind, MDR, and 
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thermal emission from the accretion disk. Wang & Robertson (1985) calculated the 
angle-integrated thermal bremsstrahlung emissivity Jb from the heated plasma in a 
propeller flow. Using their scalings, the x-ray luminosity and temperature are given by 
- AiiRl^djB ~ 4.0 X W^Bi2^l'^m%'' erg s'^ and kT ~ hQm]'^^ Bi2^'^ keV, respectively, 
where 5 is the width of the magnetospheric boundary layer where the plasma is heated by 
the magnetic fleld. This assumes spherical symmetry, so the luminosity will probably be 
~ 10'^^ ergs s~^ or less for the case of a disk geometry. If PSR 1757-24 is not a propeller 
source, we would expect the non-thermal x-ray emission which is characteristic of young 
rotation-powered pulsars ([Seward fc Wang 1988|) - e.g. a power law emission spectrum 
with a photon index F ~ 2 and an x-ray luminosity given by the empirical relation 



1 QO 

oc {VlVl) ■ ~ 1.2 X 10^^ erg s'^ At a distance of 5 kpc QGaensler fc l^rail 2000|) , the 



flux from the non-thermal, non-propeller emission would be 4 x 10~^^ ergs cm~^ s~^, which 
would be easily detectable by XMM or Chandra. The detection of dimmer, thermal x-ray 
emission instead of the brighter non-thermal emission would be evidence in support of the 
propeller model. 

Cooler disk blackbody emission could also be visible at optical and infrared wavelengths 
from the accretion disk. The spectrum of propeller disks depend on B, m, and the disk 
orientation with respect to the line of sight (e.g. [Perna, Hernquist, fc Marayan 20"U(]| ), 
but we can place an upper limit on the optical emission from a PSR B1757-24 disk in 
the following manner. Assuming that the majority of the optical disk emission originates 
at the innermost disk radii ( Perna, Hernquist, fc Narayan 2000| ), the upper limit on the 
optical luminosity is < GM^rh/ Rm ~ QAlm\^^ Bi2~'^^'^ Lq. This estimate ignores heating 
of the disk by irradiation from the pulsar, which dominates the heating only for pulsar 
luminosities > 10^^ ergs s^^ (perna, Hernquist, fc Narayan 20"00|) . Assuming B^ = 1, 



m 



1, and an apparent visual magnitude of —26.78 for the Sun (pang 1980|) , the apparent 



magnitude of the disk at 5 kpc (uncorrected for extinction) is > 20.7. The extinction 
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Ay can be estimated from the formula A'^^^ = 1.79 x lO^^Ay mag cm ^ (|PredeJil fc Schniltt| 



1995|) , where is the HI column density along the line of sight. At a distance of 5 kpc, 



Nh ~ 10^^ cm ^ in the Galactic Plane (as can be seen from the AXP spectral data in Perna. 



Hernquist, fc Narayan 20001) , which yields an extinction of Ay ~ 5.6 mag and a lower limit 
on the disk magnitude of mv > 26. This is comparable to the estimated magnitudes of 
propeller disks around AXPs as calculated by Perna et al. (2000). 

If PSR B1757-24 is indeed surrounded by an accretion disk, then fossil accretion 
disks and propeller spin-down may be present, or may have been present at one time, 
in a significant fraction of isolated pulsars. This would affect the distributions of pulsar 
ages, magnetic field strengths, and space velocities inferred using the pulsar's P and P 
values and the MDR formulae discussed in §1. Magnetic field strengths estimated using 
Equation (P, for example, would overestimate the true field strengths for both pulsars 
currently undergoing propeller spin-down and for pulsars which had experienced some 
propeller spin-down in the past (and whose disks had subsequently dissipated). The effect 
on the distribution of radio pulsar ages - and hence the distribution of pulsar velocities 
inferred from their angular positions (e.g. [Cordes fc Chernoff 1998 ) - depends on whether 



each pulsar is currently undergoing propeller spin-down or not. If a pulsar is currently 
undergoing propeller spin-down, then the MDR timing age Tmdr is an underestimate of the 
true age (as exemplified by PSR B1757-24). If a pulsar had a propeller disk which then 
dissipated, its present-day MDR age would be an overestimate of its true age, since the 
pulsar had undergone a period of increased spin-down rate (over the MDR spin-down rate) 
in the past. The latter scenario is probably more prevalent in the observed population of 
radio pulsars, since radio pulsations may be scattered and quenched in neutron stars with 
strong propeller flows (e.g. [Alpar 2000| ; [Fender et al. 19961) , and therefore one might expect 
a systematic overestimation of radio pulsar ages due to fossil disks and propeller spin-down. 
Curiously, such a systematic overestimation of pulsar ages using MDR has been inferred 
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from some pulsar population studies ( |(Jordes fc CliernoH:' 1998|) . 



4. Summary 



We have shown that the addition of torques from an accretion disk can explain the 
discrepancy between the MDR timing age of PSR B1757-24 and its true age. This model 
can be tested through x-ray and optical observations of this pulsar. The accretion disk 
model for pulsar B1757-24 leaves open the question of whether or not pulsar B1757-24 is an 
unusual and rare object or if it instead reflects a generic feature in the evolution of neutron 
stars. If the former is true, the accretion model removes the need ( Gaensler fc Frail 2000 ) 



to revise our current understanding of the physics and astrophysics of neutron stars because 
of this single pulsar. If the latter is true, however, the distributions of pulsar magnetic field 
strengths, ages, and velocities will have to be reconsidered to take into account the effects of 
increased spin-down due to accretion disk torques. Observations of young pulsars associated 
with supernova remnants may hold the key towards resolving this question, because the 
pulsar ages can be constrained independently of the pulsar temporal parameters. 
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Fig. 1. — The discrepancy between the magnetic dipole radiation (MDR) age of 16 kyr and 
the proper motion age of > 39 kyr for PSR B1757-24 may be resolved by the addition of 
"propeller" torque due to an accretion disk. This is shown by the plot of the calculated age 
(dashed lines) of pulsar B 175 7-24 versus the neutron star magnetic field strength B and mass 
infall rate m, assuming a combined spin-down torque due to both magnetic dipole radiation 
(MDR) and an accretion disk formed from supernova debris. The allowed combinations of 
B and m fall on portions of the heavy solid line (corresponding to the observed period and 
spin-down rate of the pulsar) lying outside the shaded areas excluded by upper limits on the 
pulsar proper motion (Gaensler & Frail 2000), and the condition necessary for the propeller 
effect {Rm > Re)- 



